Identification of Kelvin waves: numerical challenges 
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Kelvin waves are expected to play an essential role in the energy dissipation for quantized vortices. However, 
the identification of these helical distortions is not straightforward, especially in case of a vortex tangle. Here 
we review several numerical methods that have been used to identify Kelvin waves within the vortex filament 
model. We test their validity using several examples and estimate whether these methods are accurate enough 
to verify the correct Kelvin spectrum. We also illustrate how the correlation dimension is related to different 
Kelvin spectra and remind that the 3D energy spectrum E(k) takes the form l/k in the high-fc region, even in 
the presence of Kelvin waves. 



I. INTRODUCTION 

Dissipation of energy in the zero temperature limit is a cen- 
tral question in the field of quantum turbulence. In helium su- 
perfluids turbulence is related to motion of quantized vortices 
that generate the superfluid velocity field. At scales smaller 
than the intervortex distance energy is expected to be cas- 
caded to smaller scales via Kelvin waves (KWs) until it can 
be transformed into thermal excitations, phonons in case of 
4 He. 12 Therefore, the proper identification of Kelvin waves 
is important if we want to verify the theoretically predicted 
Kelvin spectrum 3 " 5 . Experimentally this is probably not pos- 
sible, but numerically the identification is conceivable. Ac- 
tually, several numerical simulations have tried to verify the 
correct spectrum using the vortex filament model 6 " 9 . In our 
opinion, none of those is fully convincing. The cascade to- 
wards smaller scales is evident in simulations but the correct 
value for the slope of the spectrum is dominated by a numer- 
ical noise or distorted by the method used to determine the 
Kelvin amplitude. The simulations by Kozik et al. ' , using 
the Hamiltonian formulation of the filament model, correctly 
identify the Kelvin modes, but their spectrum may depend on 
the absolute amplitude of the initial spectrum 510 . 

In BEC's the Kelvin spectrum has been determined by 
Yepez et al. 1 1 using the Gross-Pitaevskii equation. However, 
the interpretation of these results was strongly criticized 12 ' 13 . 
In high-fc region the k~ 3 spectrum (originally interpeted as 
Kelvin spectrum) most likely results from the velocity and 
density profile inside the vortex cores 13 . The low-A: region 
was originally interpreted it as Kolmogorov K41 law, even if 
the spectrum is likely visible at scales smaller than the inter- 
vortex distance (authors did not evaluate the mean vortex sep- 
aration). One option is that this spectrum comes from Kelvin 
waves, but since the fluid is highly compressible the explana- 
tions presented in the Comment 12 are more likely. 

In this article we list several ways that can be used to iden- 
tify Kelvin waves. We work in the framework of the vortex 
filament model 14 , where the vortices are considered to be thin 
and the superfluid velocity can be calculated simply from the 
vortex configuration s(| ), where £, is the length along the vor- 
tex, using the Biot-Savart integral. In other words, we con- 
sider only length scales much larger than the vortex core size, 
a. We start with two simple cases where the definition of the 
Kelvin waves is obvious: straight vortex with Kelvin waves 



and a vortex ring occupied by KWs. In both cases the Kelvin 
spectrum can be determined using a simple Fourier transfor- 
mation. We use these two sample cases to test other methods 
that have been used previously in the literature. Such meth- 
ods are, for example, curvature, energy spectrum, and fractal 
dimension. 



II. KELVIN WAVES 

Kelvin waves are helical distortions on a vortex. They 
can be generated by applying a counterflow v n — v s along a 
vortex 15 . Vortex reconnection, and more generally, interac- 
tion with other vortices, with vortex itself, or with boundaries 
also induce Kelvin waves. In the small amplitude and long 
wave length limit the dispersion relation for a Kelvin wave 
with wave vector k is given by 



Kk 2 

co k = —[H2/ka)-r\. 



(1) 



Here k = hjm\ (for 4 He) or h/2m$ (for 3 He-B) is the circula- 
tion quantum, a the vortex core size and y = 0.57721 ... is the 
Euler constant. 

On a straight vortex (taken to be along z-axis) the definition 
and identification of the Kelvin waves is simple. Assume that 
vortex can be represented as coordinates x(z) and y(z) and 
define w(z) = x(z) + iy(z). Now the vortex configuration can 
be presented as a sum of different Kelvin modes: 



w (z) =£vi^exp(ife). 

k 



(2) 



In case of periodic boundary conditions along the z-direction 
(typical for simulations), the fc-vector is discrete and given by 
k = 2nm/L z , where m (mode) is integer and L z is the period. 
For a particular vortex configuration a Kelvin spectrum can 
be defined as = \wt\ 2 + | | 2 with k > 0. Theoretical pre- 
dictions, which apply for a statistical average, predict that the 
steady-state Kelvin spectrum takes the form 3 " 5 
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with tj ps 1.5... 2. This spectrum results from interaction 
of different scale Kelvin waves and provides a constant (k- 
independent) energy flux, 0, through different scales. If e m 
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(with m > 0) determines the fraction of positive Kelvin modes, 
then we can parametrize the amplitudes using the modes: 



\w m \ = sfe^Am n , 
\w-m\ = Vl -e m Am" TJ . 



(4) 



Numerically, a great care should be used when identifying 
the spectrum. If the vortex configuration is such that the lo- 
cal curvature approaches the numerical resolution, k KS , then 
even the interpolation used to obtain the vortex configuration 
at equidistant points (needed for FFT) can result in errors that 
distorts the spectrum for £- values k > k les /5. One way to avoid 
these errors is to directly solve the Fourier coefficients. But 
then one loses the speedup obtained by FFT. An alternative 
is to keep the vortex points equally z-distributed by remov- 
ing the motion along the z-direction. By introducing a vec- 
tor p(z,t) = (x(z,t),y(z,t)), which only determines the point 
of intersection of the vortex line with the plane, one obtains 
that 1 - 16 



dp ds A ds 
~dt = dt~\ Z 'dt) l Z ' d 



dp 



(5) 



Here z is the unit vector along the z-axis. This formula, to- 
gether with the Biot-Savart equation can be used to derive the 
Hamiltonian equation for the vortex motion 116 . This projec- 
tion should be used with care, since it limits the real vortex 
motion. It, for example, prohibits a vortex to reconnect with 
itself. 

Cascade due to Kelvin waves is typically very weak. The 
energy flux, 9, can be increased by increasing the amplitudes 
of the Kelvin waves. Simulations always contain some dis- 
sipation (or noise) that typically is strongest at the smallest 
scales. For small energy flux this numerical dissipation may 
in some cases correctly mimic the energy sink. However, if 
the numerics is done in a such way that the energy is well con- 
served, this flow out from the numerical fc-region can become 
smaller than the energy flux due to Kelvin-wave cascade. This 
results in a numerical bottleneck occurs, implying accumula- 
tion of Kelvin waves near the resolution limit and generally re- 
sulting excessive fractality and noise. In simulations this can 
often be seen as a failure of the Fourier presentation for the 
vortex. Therefore, any numerical method that tries to deter- 
mine the steady-state Kelvin spectrum should contain enough 
dissipation such that this bottleneck due to finite resolution 
is avoided. How this is implemented in numerics, especially 
with the vortex filament model, is still somewhat open ques- 
tion. One option is to work at finite temperatures and have 
high enough resolution such that all the energy can be dis- 
sipated by mutual friction. Important is, that a proper dissi- 
pation is arranged dynamically. Simply damping the highest 
modes might not be enough. 

One should also note that near the resolution limit the dis- 
persion relation, Eq. (1), becomes distorted. Simple estima- 
tion states that ft) becomes flat near k res . Because the dynamics 
is not accurately modeled at the smallest scales, these scales 
can result in different spectrum and should therefore be omit- 
ted from the fits. 



III. KELVIN WAVES ON A VORTEX RING 

For a vortex ring with Kelvin waves one can also make a 
Fourier presentation. Consider a vortex ring of radius Rq in 
the z = Zo plane located symmetrically around z-axis. If we 
occupy the ring with Kelvin waves, in a similar fashion than a 
straight vortex, we may parametrize it as follows: 



x = ^R m cos(m0 + (p m ) cos 

m 

y = cos(/M0 + <Pm) sin <j) 

m 

Z = Zo- Y, R m sin(m0 + (p m ), 



(6) 



where = arctan(y/jc) is the azimuthal angle and (p m is the 
phase of the mode m. Amplitudes R m (which correspond to 
\w m \, in case of a straight vortex) can be calculated using a 
FFT in the following way: First the points along the vortex 
must be interpolated at equidistant (j) points. This can be real- 
ized using e.g. cubic interpolation. Now the Fourier transfor- 
mation of 



Z = xcos +ysin0 — iz 



(7) 



directly determines the R m 's. The zero frequency term from 
the Fourier transformation gives both the ring radius Rq (real 
part) and also the location of the ring, zo (imaginary part). 

If we want to occupy the vortex with a particular spectrum, 
we can take 



(8) 



R-m = \fl-£ m Am n . 



Here e m (with m > 0) determines the fraction of positive and 
negative Kelvin modes. 

At zero temperature, or when the normal fluid is at rest, ki- 
netic energy and momentum are only due to superfluid com- 
ponent and completely determined by the vortex configura- 
tion (we ignore externally applied velocities). Additionally, if 
vortices form loops and the vorticity disappears at the infin- 
ity, one may evaluate kinetic energy, E, and momentum, P, in 
numerically convenient way by line integrals over the vortex 
configuration: 



(9) 



E = p s K f v s -s xds 
P = ipj? 



(10) 



These formulas make it possible to evaluate the accuracy of 
the numerical scheme used. At zero temperature both the en- 
ergy and momentum should be conserved. However, even if 
the energy is well conserved, e.g. by accuracy of 0.1 percent, 
this fluctuation is still very large if we compare it with the en- 
ergy related to smallest scale Kelvin waves. Therefore, the 
errors coming from the time iteration scheme can be large at 
the smallest length scales. In the following we are not concen- 
trating on these errors but pay attention only to methods that 
try to identify the Kelvin waves on a tangle. In other words 
we assume that the vortex configuration is correct at the nu- 
merical points s,-, i = 1, . . . ,7V. 
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ously observed (with somewhat limited resolution) to take a 
form c~ a . 9 At least our test with spectrum given by Eq. (2), 
indicate that this kind of fit is poor (exponential fit is often 
better) and the obtained a is very sensitive to the fit region. 
More importantly, the location of the maximum and also the 
shape of the histogram at high curvature regions depends on 
the resolution for spectra tj < 2.5, as seen in Fig. 1. 

In case of straight vortex with small amplitude Kelvin 
waves one can relate the curvature spectrum with the Kelvin 
spectrum in the following way: 3 



\ C ^)\ 2 d^ (\w"(z)\ 2 )d Z : 



E* 4 M 2 . 



(11) 




FIG. 1. (Color Online) Curvature histogram obtained for a straight 
vortex with Kelvin waves. Histograms are averages over 50 different 
sets, where the phases and the fraction of positive and negative modes 
is randomly set. Top: Curvature histogram for a Kelvin spectrum 
with 77 = 1.7 and having amplitudes A/L z = 0.01, 0.05, 0.1, 0.25. 
Number of points is N = 8192 and all the modes with \m\ > 5 are 
present. Bottom: Curvature histogram for A = 0.1L ; . The solid (red) 
line is calculated with N = 8192 points and using all modes (up to 
the resolution). The dash-dotted (black) line is calculated with four 
times better resolution (N = 32768), also using all modes. The dashed 
(blue) line is also with N = 32768, but having only the same modes 
than in the low resolution run. The main figure is for 77 = 1.7 and the 
inset is for 77 = 2.5. 



IV. CURVATURE 

In case of vortex tangle the identification of Kelvin waves 
becomes tedious. One easily seeks different indirect meth- 
ods to identify small scale Kelvin waves. The most obvious 
thing is to calculate local curvature c(£, ) = |s"| along the vor- 
tex. An increasing average curvature is a clear sign that more 
small scale structures appear. A histogram describing the dis- 
tribution of curvature gives some more information, but our 
tests indicate that curvature histogram is a poor method to de- 
termine the Kelvin spectrum (exponent 77). For example, the 
location of the histogram maximum depends, not only on the 
spectrum, but also on the amplitude, as illustrated in Fig. 1. 
The location of this maximum, as a function of amplitude, 
is also non-monotonous. Similar effect is also seen if one 
changes tj. 

At high curvature regions the histogram has been previ- 



Therefore one might expect that the curvature spectrum could 
reveal the Kelvin spectrum, also in more complicated tangles. 
We observed that the curvature spectrum follows this law only 
when the amplitude is small and 77 > 2, i.e. only when the 
curvature spectrum °= k y is such that y < 0. For any smaller 77 
the curvature spectrum becomes flat and nothing can be said 
about 77. Improving the numerical scheme for determining 
the local curvature does not resolve this problem. If the Kelvin 
amplitude is not small, the above relation, Eq. (11), is satisfied 
even less accurately. 



V. ENERGY SPECTRUM 

Energy (or velocity) spectrum E(k) describes the distribu- 
tion of the kinetic energy at various length scales. At absolute 
zero temperature the total kinetic energy is given by super- 
fluid component^ = Ek, s = hps! v 2 d 3 r = J E(k)dk. Where 
now k is the three-dimensional vector in the momentum space 
and k = |k|. (Note that in case of Kelvin waves the ^-vector 
was ID object.) Assuming a quantized vortex with singular 
distribution of vorticity and that the vorticity disappears at the 
infinity one may derive following formula for the kinetic en- 
ergy caused by the vortices 17 : 



Ets — 



2(2n) 



d\ 



^ 2 eik-(s,- S 2) (g ;. § ^ (12 ) 



Here the double integration is along vortices, described by co- 
ordinates and ^2- Notations s,- = s(|,-) and §■ = c/s(<^)/c/<^ 
are the location of the vortex core and tangent at respec- 
tively. Additionally one should emphasize that this formula- 
tion is exact only if the vortices form closed loops and that 
the vorticity disappears far away. Otherwise, at least the low 
fc-values are calculated incorrectly, which is not emphasized 
inRefs. [17-19]. 

Fortunately, the integration over different ^-directions in 
Eq. (12) can be done analytically 20 , resulting that 



E{k) 



PsK~ 
{2nf 



d& / di; 



sin(fc|si -s 2 | 
1 fc|si-s 2 | 



■(§'1 



§2) 



(13) 



For anisotropic situations, the dependence on the absolute 
value of the wave number, k, should be understood as an angle 
average. 
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In our numerical scheme vortices are described by sequence 
of points and we assume that between these points vortex is 
straight. This is a standard assumption used also by several 
other authors. The numerical integration of the above double 
integral is quite straightforward but, due to the oscillating part 
of the integrand, some care should be taken when k is large. 
However, we found that it is important to check this limit, be- 
cause the energy spectrum at fc-values that approach and ex- 
ceed the numerical resolution limit should converge towards 
k^ 1 spectrum. This is not the case in some of the previous 
vortex filament calculations presented so far 17 ' 20 ' 21 and there- 
fore some caution should be kept in mind when reading those 
papers. Additionally, one numerical check for the spectrum is 
given by the absolute amplitude of the spectrum: in the high-fc 
region, near the resolution, the amplitude should be given by 
the vortex length. 

For a single vortex ring with radius R, the line integration 
formula, Eq. (13), can be simplified further. Simply putting 
Si = R(cos ; x + sin0,-y) when §'j • s' 2 = cos(0i — 02 ) one ob- 
tains, after few standard integration techniques, that 

p,K 2 R 2 f K , sm(2kRsmx) 

E(k) = — / dx — i ^cos2x . (14) 

w K J 2kRsinx y ' 

This can be presented using generalized hypergeometric 
functions, but we found that it is much more convenient to 
do the integration numerically. The spectrum for vortex ring 
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FIG. 2. (Color Online) Top: Energy spectrum E(k) of vortex ring of 
radius R. At large scales (kR <S 1) E(k) = p s K 2 R A k 2 /d and at small 
scales (kR > 1) E(k) = p s K 2 R/(2k), indicated by the dashed lines. 
Bottom: Deviation of the energy spectrum from the kR 2> 1 limit. 
The amplitude of the oscillations scales like k~^l 2 . 




FIG. 3. (Color Online) Energy spectrum E(k) of a straight vortex 
with Kelvin waves. The spectrum for Kelvin waves is given by rj = 
1.7 with amplitude A = 0.1L-, where L z = 1 mm is the period along 
z-direction. The two straight lines corresponds spectrum of straight 
vortex with length L = L- (lower) and L = 1 .4347L- (upper), which 
is the total length of the vortex (per period used). 



is presented in Fig. 2. At large scales the energy spectrum is 
given by E{k) = p s K 2 R 4 k 2 /6 which can be obtained by doing 
a series expansion at small k for Eq. (14). At small scales the 
vortex ring looks like a straight vortex (of length 2nR) and the 
energy spectrum is given by E{k) = p s K 2 R/(2k). The spec- 
trum for single ring is very similar to the one produced by four 
rings in Ref. [21], but our better resolution makes it possible 
to observe small amplitude oscillations. We found that these 
oscillations, with period n/R, around this limiting curve scale 
like A:~ 3 / 2 , which is illustrated in Fig. 2. To resolve these os- 
cillations the discretization step for different fc-values should 
be smaller than this oscillation period, which makes the cal- 
culations of high fc-values time consuming. 

An advantage of the line integration formula over the stan- 
dard Fourier transformation used e.g. by Kivotides et al? x is 
that one avoids the numerical problems that is caused by the 
divergent 1 jr velocity field generated by vortices. A disad- 
vantage of this method is that in the presence of the walls or 
with periodic boundaries some care should be taken since not 
all of the assumptions in deriving Eq. (12) are fulfilled. Typ- 
ically the boundaries result in extra surface terms that make 
the correct calculation more cumbersome. For a single plane 
boundary the energy spectrum can be calculated by extending 
the second line integral in Eq. (13) to include also the image 
vortices. Nevertheless, the high fc-values are correctly taken 
into account even with this simple formulation. 

Even if the 3D energy spectrum E(k) is suitable for deter- 
mining the large scale velocity structures appearing in vortex 
tangle, it is not sensitive for small scale Kelvin waves that 
have characteristic scale smaller than the intervortex separa- 
tion. This is illustrated in Fig. 3 where we have plotted the 3D 
energy spectrum for a straight vortex with and without Kelvin 
waves. The l/k contribution coming from the straight vortex 
is so dominant that very little can be said about the Kelvin 
spectrum, at least for spectra with 77 > 1.5. The increased 
length only raises the spectrum up without affecting the slope 
in the high-£ region. In case of vortex ring the identification 
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of Kelvin waves is even more complicated due to the oscilla- 
tions that are already present for a ring without Kelvin waves. 
This just emphasizes that in order to identify the Kelvin spec- 
trum the direct measurement of the vortex position becomes 
more important than determination of the 3D energy spectrum 
at high fc-values. 



VI. FRACTAL DIMENSION 

At low temperature vortices typically looks quite wiggly, 
almost fractal like. This is due to Kelvin waves. Therefore, 
one way to characterize the tangle is to determine the fractal 
dimension for vortex lines 9 . 

Upper-bound estimate for the fractal dimension can be cal- 
culated by determining the correlation dimension D intro- 
duced by Grassberger and Procaccia 22 . In case of N dis- 
cretization points, D can be determined by calculating the 
number of point pairs, «, whose separation is less than e. In 
the limit of e — > the correlation integral K(e) = n/N 2 takes 
the form^e) e D . 

How are different Kelvin spectra and D related? In order to 
answer this we have determined D in case of vortex ring that 
is occupied by Kelvin waves with different spectra and with 
different amplitudes. Figure 4 summarizes our results. As 
expected, D is clearly bigger than unity only when rj < 1 .5 
and when the amplitude is large enough. This limits the us- 
ability of the fractal dimension to spectra that are less steep 
than predicted theoretically. The dependence on the ampli- 
tude makes it difficult to determine 77 accurately. However, 
if one observes a correlation dimension that is clearly larger 
than unity, then one can almost certainly state that 77 < 1.5. 
Also, even if we used a simple vortex ring to determine the 
fractal dimension, the results should be valid more generally 
at scales smaller than the average intervortex separation, i.e. 
in the region where the Kelvin cascade is important. 

The difficulty in determining the fractal dimension is that 
numerical noise near the resolution limit (appearing because 
there are too few pairs) causes that D — » 3. We observed that 
this happens especially when the amplitude is large and 77 is 
small (near unity). Therefore, the fit region should be limited 
to scales larger than the numerical resolution. The incorrectly 
chosen fit region might explain why the calculated correlation 
dimensions (of order 1.5) in Ref. [9] are much bigger than the 
presented configuration (Fig. 20 in that article) would suggest. 
For example, Kelvin spectrum with 77 = 1 .20 and with ampli- 
tude A = Q.25R, presented in the inset of Fig. 4, results in a 
fractal dimension of D ps 1 .20 only. 



VII. SMOOTHED VORTEX 

For a vortex tangle the identification of Kelvin waves can by 
done by using the concept of "smoothed vortex" 1 ' 9 . From the 
original discretized vortex filaments, determined by points s, 
(i = 1, . . . , AO, one can construct a smoothed vortex line s smoot h 
by using every nth points as nodes for a cubic-spline interpo- 
lation. One can then define the Kelvin-wave amplitude a(^) 




A/R 



FIG. 4. (Color online) Correlation dimension D as a function a 
Kelvin-wave amplitude for a vortex ring occupied by Kelvin waves 
\m\ > 6 for different spectra |w m | <* with 77 = 1.00, 1.05, ... 1.50 
and using 4096 points. Calculated D is an average over several sam- 
ples where the phase and the distribution of positive and negative 
Kelvin modes is randomly set. The inset shows one configuration 
when 77 = 1.2 and A/R = 0.25. 



as the distance between the smoothed and the original fila- 
ment: a(^) = \s — s smoot h|. The amplitude spectrum A (A:) can 
be defined as 

l -Ja 2 {^=J^A{k)dk. (15) 

We have applied this method for a straight vortex with Kelvin 
waves, where the spectrum takes some known form. Figure 5 
illustrates that the above scheme has a tendency to follow the 
correct spectrum. For type A (k) <^ k~P spectrum, /3 obtained 
from the fit is typically too big for spectra with j3 < 4, and too 
small for spectra with /3 > 4. In some cases the method can 
result in a spectrum that is totally incorrect. This is also shown 
in Fig. 5. Better identification of Kelvin waves could perhaps 
be obtained if the smoothed vortex was allowed to avoid the 
original data points. This removes the discontinuities of the 
derivatives a'(^) appearing when a(E,) = 0, which is one rea- 
son why the spectrum tends towards /3 = 4. We have tried 
this type of interpolation but so far the improvements are only 
minor. 



VIII. CONCLUSIONS 

Here we have presented several methods that can be used to 
identify Kelvin waves. In case of a straight vortex or a vortex 
ring the identification can be done most reliable. However, in 
case of vortex tangle a great care should be used. The biggest 
problem is the lack of proper definition of a Kelvin wave on a 
curved vortex. In order to properly define a Kelvin wave there 
must exist a scale separation between the characteristic size of 
the underlying vortex and its disturbances. 
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FIG. 5. (Color online) Kelvin-amplitude spectrum A(k) calculated 
using a "smoothed vortex" -method described in the text, with 1024 
points and n= 16. The (red) solid lines describe the correct spectrum 
and the dashed lines are fits at large k, where the method should be 
applicable. Top left: spectrum with = 2 (slope for the fit is 2.4). 
Top right: spectrum with [3 = 3.4 (fit 3.5). Bottom left: spectrum 
with j8 = 5 (fit 4.0). Bottom right: spectrum with /3 = 3, but omitting 
all modes with m > 100 (fit gives 4.1). 



Generally the Kelvin cascade to smaller scales can be qual- 
itatively identified using, e.g. average curvature. An accurate 
determination of the Kelvin amplitude is still very demand- 
ing and can result in large errors, depending on the numerical 



method used. It is alarming that some numerical methods can 
result in a spectrum that is totally incorrect but still very close 
to the spectrum predicted theoretically. Therefore, all methods 
should be properly tested using a simple configuration where 
the spectrum is known in advance. 

The highly dominating geometrical 1 /k scaling in the 3D 
energy spectrum at high-fc illustrates that at small scales the 
determination of the vortex location becomes much more es- 
sential than the knowledge of the energy spectrum. Only af- 
ter determining the vortex location, we have some hope for 
identifying the Kelvin spectrum. Identifying the vortex loca- 
tion is not a problem with the vortex filament model and is 
quite straight forward with the Gross-Pitaevskii equation. Ex- 
perimentally the identification is challenging but recently the 
visualization of quantized vortices has become possible 23 . 

Furthermore, the accuracy of the numerical schemes typi- 
cally limits the proper determination of the Kelvin spectrum 
to scales k < k KS /5, or even less. Similar requirement results 
from the dissipation, which must be present in order to 
avoid numerical bottleneck that would otherwise generate 
in a non-physical fractalization of the vortex at the smallest 
scales. Therefore, the identification of Kelvin waves on a 
tangle is a great challenge and requires resolution that is 
much higher than the average vortex separation. Here faster 
computers and especially new numerical algorithms, like the 
tree method 24 , are of great importance. 
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